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of the entanglement of the initial coin state. For states with a special entanglement 
structure it is shown that this entanglement can be meausured with the mean value 
of the walk, which depends on the i-concurrence of the initial coin state. Further on 
the entanglement evolution is investigated and it is shown that the symmetry of the 
probability distribution is reflected by the symmetry of the entanglement distribution. 
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1. Introduction 



Quantum walks are the quantum counterparts of classical random walks. Since they 
were introduced in pp, they were investigated to construct the same powerful algorithms 
for future quantum computers as in the classical case. An overview can be found in [2j. 
Here we are interested in an entanglement measure and will study a special constructed 
discrete quantum walk scheme with more than one coin. We will connect the 
entanglement of the initial coin state with the outcoming mean value in position space. 
An equation will be given that establishes the mean value as an entanglement measure 
for the initial coin state. Further we will investigate the evolution of the entanglement 
in coin space. 



1.1. Quantum walk scheme 



We use a simplified quantum walk scheme with M coins, which is in it's base similar to 
the models used by Flitney et al. [3j and Brun et al. p!. 

The Hilbert space consists of two parts, a (M x 2)— dimensional coin space and a 

n— dimensional position space Hp. The quantum walk scheme consists only of applying 
repeatedly a special transformation operator E to an inital state = 0)) = |^o)- 

|^(t)) = E*|^o) (1) 

The wave function |\l/o) consists of two parts, a part describing the coin state and 
the other giving the probability of the walker to be at a certain lattice site, |\l/o) = 
\ipc) ® \ipp)- The transformation operator E has two parts, one for the coin tossing 
and one for the shifting of the walker on the lattice. As an example we show here an 
operator with M = 3 coins, which executes the operation depending on the state of 
qubit 2: 



E 



-1 



V 



V qubit 2 



\ 



(2) 



U 



(3) 



qubit 2 qubit 2 

We use the unbiased tossing operator 

which means that after one step we have the same probability for going one step to the 
left or one step to the right. The shifting operators move the quantum walker on the 
one-dimensional lattice, = \x + and S~ = \x — l){x\, where x marks 
the position. The I{p,c} are unity operators in the postition resp. the coin subspace. 
One can also apply the operation to qubit 1 resp. qubit 3. By introducing more Ic 
operators in coin space, one can extend the E operator to more coins. 
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2. Analytical evaluation 

2.1. Fourier transformation 

We use the methods described by Nayak and Vishwanath jl] and by Brun and coworkers 
[Sj to evaluate our scheme analytically. After applying a discrete fourier transformation 
to the quantum walk scheme, we get the U operator in k-space: 




(4) 



The time evolution can then easily be calculated by the eigenvalue decomposition of the 
Uk operator. The calculation of the eigenvalues yields Ai,2 = :^(cosA; ±iVl + sin^ k). 
If we apply a similarity transformation with 

(5) 




and c± = e*'^(sin + sin^ k) we can easily calculate the matrix product of the Uk 

operator: 

with the abbreviations a = cost9 + i ^ sint6', b = , sint6', Q = arccos 

y l+sin^ k y l+sin^ k V i 

with |ap + |6p = 1. 



2.2. Calculation of the moments 



Brun et al. jH] derived an equation for the calculation of the moments of the distribution 
in position space 



dk (MUD 



dk"^ 



Ut 



(7) 



which we will use to calculate the mean value (m = 1) and the variance (m = 2). \(j)o) 
is the initial state in coin space. We start with the simplification of the integrand for 
the calculation of the first moment: 



a' b' 

{-b*y {a*y 



a'a* + b{b*y a*b' - {a*yb 
a'b* - a{b*y a{a*y + b'b* 

We find a special property of Ci since 

Ci + cl = a* a + bb'* + a*a + b'b 

= iaa*y + ibb*y=i\a\^ + \b\'') 




(8) 



2\/ 



(9) 
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Figure 1. The numerical results of the integral ci as a function of time and linear fit 
with ci{t) = aot + fli with oq — —0.2932 and ai ~ 0.5116. 



and therefore ^ J^T^dk (ci — c^) = 0. In the following we will use the abbreviation 

ci:=:^ [ dk ci. (10) 

From the symmetry of the integrand we can further conlude that ^ J^^ dk (rfi — c/*) = 0. 
Our results can easily be verified with these two simplifications. In Fig. [T]we evaluated 
numerically the integral as a function of time and made a linear regression of the 
result. 

For the calculation of the variance and of the mean deviation one needs the second 
moment. For this we evaluate the integrand of ((Tj) for m = 2. 

^ = f t 1 I . I = (11) 




a"a* + b{b*)" a*b" - 
a"h*-a{h*)" a{a*)" + b"b* 




(12) 



We know from numerical evaluation, that the integral over the matrix {UlY Ul is 
diagonal and that the entries have the same value: — ^ J^^ dk (c2 — C2) = The integrals 
over the nondiagonal entries vanish. With these two informations we can conclude that 
the second moment does not depend on the initial state, but only on the integral 

C2 := -— / dk C2. (13) 

This is in accordance with the results obtained by Konno [HIIZ!- In Fig. El we evaluated 
numerically the integral £2 as a function of time and fitted the result with a quadratic 
dependence. 

We now apply our quantum walk scheme to five initial coin states, which differ in 
the entanglement structure. It is shown that for states with a certain entanglement 
structure, the following equation for the squared mean value is vahd: 

{xr,=c^\t){l-IC^) (14) 
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Figure 2. The numerical results of the integral C2 as a function of time and fit with 
<5(i) = bot'^ with bo = 0.35i249. 



Table 1. Entanglement structure of the example states. 2-qubit entanglement as 
measured by the concurrence I1U| . 3-qubit-entanglement for tripartite states as 
measured by the tangle and 3-qubit resp. 4-qubit entanglement for fourpartite 
states as described in With pure entanglement we mean, that the state contains 
only one kind of entanglement. 



state 


eq. 


2-qubit ent. 


3-qubit ent. 


4-qubit ent. 


pure ent. 




\jGHZ) 




no 


yes 




yes 


yes 




m 


yes 


no 




yes 


yes 


i\^r) + \^8))/V2 


m 


yes 


yes 




no 


no 




(El 


yes 


no 


no 


yes 


yes 




(EH) 


no 


yes? 


yes? 


no 


no 



{x)i is the mean value related to the E operator which acts on qubit i. ICi is the i- 
concurrence [S] related to the i-th qubit of the initial coin state. In addition the following 
equation for the variance can be derived: 

{{x') - {x)')^ = (cl(t) - c,\t)) + c,\t)ICf (15) 

The i-concurrence measures the entanglement between two subsystems A and B and 
can be written as ICa-b = a/2[1 — Tr{p\)], with the reduced density matrix pA- We 
use the notation ICa-b = ICa- Our main equation (fT^ connects the mean value of the 
quantum walk with this entanglement measure for the initial state. We show analytical 
and numerical examples that eq. (fT^ only holds for so called pure entangled states, 
that means states with only one kind of entanglement, e.g. for a tripartite state either 
2-qubit entanglement or 3-qubit entanglement. For states with mixed entanglement the 
mean value stays zero. 

3. Application to example states 

We will apply our scheme to pure 3- and 4-qubit initial states, which differ in their 
entanglement structure. This known entanglement structure of these states is described 
in [12] and summarized in Table [T] We start with a parameter dependent GHZ jTH] 
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Table 2. Parameter Aq of the fit of the squared mean values, with {x)f ~ Ao{l — ICf), 
for t = 50. |0i)i resp. \(f>i)2 means the state dependent on one parameter, when 
the other one is constant. 



state 


ICi 


IC2 


IC3 


Id 


\-fGHZ) 
\A) 


202.634 
202.631 


202.634 
202.631 


202.634 
202.631 




l<^i)i 


202.639 
202.641 


202.636 
202.635 


202.639 
202.641 


202.634 
202.633 



state 



I7GHZ) = 7|000) + v^T^^llll) (16) 
which is genuine tripartite entangled for 7 g]0, 1[. The i-concurrence is the same for all 

^2 

'{1,2,3} 



reduced qubits, /C|^2 3} — ^7^(1 — 7^)- The analytical solution for the mean value is 



found to be 

(^){l,2,3} = c~i^(27^ - 1)^ = ci\l - ICl2,3}) (17) 
and we have proven the result given in (|14j). In Table |21 we give additionally the fit 
parameter Aq which comes out of the fit of the simulated mean value with {x)f = 
Aq{1 — ICf). As one can see the simluation and the analytical value fit very well. 
The state l^pe) is a pure 2-qubit entangled parameter dependent eigenstate of a certain 
spin chain and has the following form: 

l^e) = 'tilOOl) + KslOlO) + KillOO) (18) 

with the norm 2k1 + = 1 and the abbreviations Ki = and K2 = — with 

rj := a/12 + A(A — 4) and x ''? + ^^2. The i-concurrence can be calculated in terms 
of the parameters, IC^i = 1 — ^2 resp. /C| = 4(/t2 — ^2)- Again it can be shown that 
the mean values fulfill eq. ()14j) : 

(^){i,3} = 4 = ci' (1 - ICl^3}) (19) 

{x)l =c,'{l-24y = c,'{l-IC^) (20) 

In Fig. iniwe show the parameter dependence of the squared mean value and the result 
of the fit. 

The next state at which we will have a closer look is a superposition of two states: 

H^y) + \i^s))/V2 = fKilOOl) + /€2|010) + /cilOll) + (21) 



Ki|100) + K2IIOI) + fi;i|110)j/V2 (22) 

with abbreviations and norm from above. As we have shown in P2J; the entanglement 
structure of this state shows 2-qubit as well as 3-qubit entanglement. The i-concurrences 
for this state could easily be calculated, but are not so important in this case, because 
the mean value is 0, independent from any parameter: 

(a;){i,3} = {ci + cl) + 2KiK2(rfi -dl) = 

{x)2 = {ci + cl) + 2K\{di - d*) = (23) 
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Figure 3. Parameter dependence of the squared mean value for the state \tpQ) for 
time t = 50. The dots mark qubit 1 resp. 3, the x mark qubit 2. The hues are fits 
with {x)l ^ Aq{1~ ICf). 

These calculations are consistent with numerical simulations. Thus it is confirmed that 
for mixed entanglements the mean value stays at zero. 

The two 4-qubit states |0i) and |02) are also eigenstates of a spin chain [12]. The 
parameters ctj and depend on two further parameters. The state |0i) 

|0i) = aillllO) + aallOll) + aglOlll) - aslllOl) (24) 

with the norm a\ + a\ + 2a| = 1, is only twopartite entangled. For simplicity we will 
only regard the i-concurrence reduced on qubit 1 and 3, JC|^ = 4(^3 — o.'^)- It is again 
nicely seen that the mean value can be described with our equation ()14|1 : 

(^){i,3} = (2al - l) ' = cl^ (l - /q\3}) (25) 

For the two other mean values can also be shown that they are connected to the i- 

concurrences, via eq. lfT^ . 

The entanglement structure for the state |02) 

102) = -/3i|0011) + AlOllO) - AllOOl) + /3i|1100) - /^slOlOl) + /?2|1010) (26) 

with the norm 4/?^ + 2/5| = 1 is a more complex state, since there is no possibility to 
show genuine 3-qubit or 4-qubit entanglement in a fourpartite state. But in ^2] we have 
shown that this state has besides the 2-qubit entanglement either 3-qubit or 4-qubit 
entanglement. The i-concurrences reduced to one qubit are parameter independent, 
IC^i 2 3 4} = 1- It can be shown for this state, that the mean values are 0, 

(x){i,2,3,4} = {Wl + f^l) (c~i + cl) = (27) 

These results are consistent with the 3-qubit problem and therefore we suggest the 
following general result, although this is not a general proof: For purely entangled 
coin-states (2-qubit, 3-qubit,...) we expect a square mean value proportional to the 
corresponding i-concurrence, while for mixed entanglements the squared mean value is 
always zero. 
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Figure 4. Global entanglement Q as a function of time for the 7GHZ state as initial 
state, with 7 = 0.3, viewed from different lattice points. The straight line is for x = 60, 
the dashed line for x = 50 and the dotted line for x = 40. The starting point is at 
X = 60. 



4. State evolution in coin subspace 

To clear this point further we look in the following at the entanglement evolution of the 
wavefunction in the coin subspace. The complete wave function at time t consists 
of two parts, one part in position space and additionally a part in coin subspace, 

= J2x 1-^) ® l'^Coin('^))- ^^^^ have a closer look at the states |V'Coin(-^))- 
With an artificial normalization of these states at each site one can calculate the global 
entanglement measure Q |14j . 

We note the following results from our simulations: a) For all starting coin states the 
entanglement structure is constant in time at the starting point on the lattice; b) at 
other lattice points we observe entanglement oscillations. These two effects are shown 
in FiglH with the 7GHZ state as initial state. If we sit on the starting point at a; = 60 
the entanglement as measured by the global entanglement Q is constant in time. If we 
go to X = 50 or X = 40 we can observe the described entanglement oscillations. 

Another effect is shown in Fig. El The entanglement distribution over the lattice is 
plotted for t = 50 and the 7GHZ state as initial state for two different values of 7. For 
7 = 1/ the mean value in position space is and the entanglement distribution over 
the lattice is symmetric. If we take 7 = 0.3, the mean value is not equal to 0, and the 
entanglement distribution is asymmetric. 

5. Conclusions and outlook 

In conlusion, we have shown that for initial states with a special entanglement structure, 
the resulting mean value of a multiple coin quantum walk scheme measures the entan- 
glement of the initial coin state. We assume that our proposed connection is valid for 
what we call pure entangled states, states with only one kind of entanglement. Further 
on it is shown, that the symmetry of the probability distribution is reflected by the 
entanglement distribution. 



Entanglement measurement with discrete multiple coin quantum walks 



9 



O" 0.5 




site 



100 120 



Figure 5. Distribution of the global entanglement Q over the lattice for t — 50 for 
two different 7— paramters for the 7GHZ state as initial state with 7 l/\/2 (dashed 
line) and 7 — 0.3 (straight line). 
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